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Abstract 

The numerical results for the computed moduli of the irreducible three-loop con- 
tributions to the thermodynamical pressure of an SU(2) Yang-Mills theory in the 
effective theory for the deconfining phase are explained in detail. Irreducible three- 
loop integrations are compared with two-loop integrations and the different nature of 
their integrations is scrutinized and illustrated numerically. The numerical results 
show a rapid convergence in the loop expansion of Yang-Mills thermodynamics. 
The statistical method used for irreducible three-loop integrations is explained and 
checked for two-loop integrations. The statistical results for two-loop integrations 
are compatible with the former computed analytical results showing the reliability 
of the statistical method. This is a companion paper to [1] . 



1 Introduction 



One of the major problems in constructing Yang-Mills thermodynamics is 
that a reliable approximation of high-temperature thermodynamics related 
to 4D Yang-Mills theories in terms of the small-coupling expansion seems to 
be impossible. 

The nonconvergence of the small-coupling expansion is due to the fact 
that an empty (trivial) ground-state is invoked to construct an approximat- 
ing series for the full partition function. Meanwhile, fluctuations of nontrivial 
topology have a profound impact on the ground-state estimate, and are com- 
pletely ignored in small-coupling expansion since their weight posesses an 
essential zero at vanishing coupling. 

Consequently, the strong correlating effects of these extended field config- 
uratons are completely ignored. This is a fact which is expressed by tree-level 
masslessness and only week radiative screenings of all gauge bosons leading 
to the nonconvergence of the expansion. However, by considering an a priori 
estimate for the ground state of an SU(2) Yang-Mills theory at high tem- 
peratures, which is obtained by a self-consistent and sufficiantly local spatial 
coarse graining over interacting and stable BPS saturated topological field 
configurations, according to the effective theory reviewed in [1], leads to a 
rapidly converging loop expansion [H [2] . 

The argument of [2] is that a self-consistent spatial coarse-graining, which 
involves interacting (anti)calorons of unite topological charge moulus, im- 
plies that real time loop expansions of thermodynamical quantities in the 
deconfining phase of SU(2) and SU(3) Yang-Mills thermodynamics are, IPI 
resummations, determined by a finite number of connected bubble diagrams. 

To see this in more detail, recall the following fundamental aspects of 
the effective theory. In the effective theory a composite adjoint Higgs field 
which describes the topologically nontrivial part of the ground state is in- 
troduced and implimented. This field is associated with periodic instantons 
(Calorons of topological charge one and trvial holonomy) of high temperature 
in the deconfining phase and is used as a background for coarse-grained topo- 
logically nontrivial sector of the theory. This field is quantum mechanically 
and thermodynamically stabilized. 

A macroscopic pure-gauge configuration which is the solution of the equa- 
tion of motion for the topologically trivial sector in the presence of the back- 
ground is implimented in order to include the interactions between trivial 
holonomy Calorons. 

As the modulus of the Higgs field decreases with the temperature as 

101 ~ w where A is the Yang-Mills scale, the effects of topological defects 
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die off at large temperature in a power-like fashion. Asymtotic freedom 
and Infrared-ultraviolet decoupling of the fundamental thery, which are the 
results obtained in perturbation theory at zero temperature, are preserved. 

Thermo dynamical quantities are in this framework calculated as loop 
expansions about the nontrivial ground state consisting of the Higgs field and 
the pure gauge configuration. The tree-level excitations are either massive 
thermal quasiparticles or massless photons, which interact very weakly. The 
effective theory is both infrared-and ultaviolet finite. The former property 
is due to the existence of Caloron- induced gauge boson masses (IR cut-off), 
and the later to the compositeness constraints on the loop momenta arising 
from the existense of the composite scale of the Higgs field (UV cut-off). 

By taking the topologically nontrivial contributions through the effective 
theory for deconfining SU(2) Yang-Mills thermodynamics into account, the 
composite Higgs field induces a dynamical gauge symmetry breaking from 
SU(2) — > U(l) implying that two out of the three propagating and coarse- 
grained gauge modes acquire temperature dependent mass. According to this 
and the fact that the off-shellness of these modes, along with the momentum 
transfer in local vertices, are highly constrained by spetial coarse-graining a 
rapid convergence the loop expansion is expected [H [2] . 

The aim of this paper is to provide numrical evidence for the conver- 
gence in the loop expansion of Yang-Mills thermodynamics [Tj. It is worth 
emphasising that this paper contains illuminating details and some technical 
aspects not present in pQ. 

The paper is organized as follows. Section 2 considers the effective gauge 
coupling, which lies at the heart of the numerical results, and the relating 
Feynman-rules with the constraints on the momenta for loop integrations 
are also explained. These were not displayed explicitly in [T]. Section [S] 
compares irreducible three-loop diagrams with two-loop diagrams, arguing 
analytically that integrations concerning the former are, in contrast to the 
later, either compact or empty and therefore give much more suppressed con- 
tributions to the thermodynamical pressure. Section 4 gives a mini review 
of the Monte-Carlo method used for loop integrations and verifies numer- 
ically the analytical claim of section 3 by using numerical illustration and 
Mote-Carlo integration. Numerical illustration shows that the region of ra- 
dial loop integration for irreducible three-loop contributions is, in contrast 
with the region of radial loop integration for two-loop, either compact or 
empty. Monte-Carlo integration shows that irreducible three-loop integra- 
tions lead to much more suppressed contributions to the thermodynamical 
pressure than two-loop integrations. This is precisely in agreement with the 
result of the analytical argument of section 3. The stability and reliability 
of the Monte-Carlo method used for doing loop integrations is also tested by 
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comparing the results of the Monte-Carlo integration for two-loop with the 
former analytical result for two-loop. Section 5 summerizes the objectives 
achieved and concludes with an open problem. 
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2 Effective gauge coupling and Feynman rules 



In [T] the effective theory was reviewed by taking the effective action as the 
startig point resulting into an evolution equation A(a) for temperature as a 
function of tree-level gauge boson mass. This evolution has two fixed points 
a = and a = oo, where the lowest and highest temperatures Ac = A(a = oo) 
and Ap = A(a = 0) are attainable, respectively. 

The evolution A (a) can be inverted to yield an evolution e(A) for the 
effective gauge coupling as a function of temperature. In what follows the 
evolution of e(A) is displayed, since it lies at the heart of the numerical 
results, and the relating Feynman rules, which depend on the effective gauge 
coupling, are also written down explicitly. These were used in [T] for doing 
loop integrations, but were not presented explicitly. For SU(2) and SU(3) 
this is illustrated in figure 2.1. The evolution of e with temperature exhibits 
a logarithmic pole, e oc — log(A — Ac), where Ac = 13.89 denotes the critical 
value of the dimensionless temperature A = and the value of e at the 
plateau is e = v^vr ~ 8.89.. 



e 




2.5 : 

100 200 300 400 500 ^ 

Figure 1: Evolution of the effective gauge coupling e in the electric (decon- 
fining) phase for SU(2) (grey line) and SU(3) (black line). 

In total, the theory seems to have three phases: The electric phase at 
high temperatures, the magnetic phase for a small range of temperatures 
comparable to the Yang-Mills-scale A and a center phase for low tempera- 
tures. The electric phase is deconfining, the magnetic phase is preconfining 
and the center phase completely confining. Here we are only interested in 
the deconfining (electric) phase. 
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Now by knowing the effective gauge coupling one can formulate the Feyn- 
man rules in the unitary- Coulomb gauge using the real-time formulation 
of finite-temperature field theory. The real-time is preferable because the 
implementation of constraints on the momenta (([6]), (171), ([9])) is rather in- 
convenient in the imaginary-time formalism. It should be noticed that the 
unitary- Coulomb gauge is a completely fixed gauge, therefore no Faddeev- 
Popov determinants need to be considered and no ghost fields need to be 
introduced. According to the discussion in the introduction the effective the- 
ory has a stabilized, composite and adjoint Higgs field (f) characterizing its 
ground state, where in the unitary gauge (f) is diagonal and the pure-gauge 
background is zero. 

The physical gauge choice for a residual gauge freedom due to the unbro- 
ken Abelian subgroup U{1) is the Coulomb gauge. In the unitary-Coulomb 
gauge each of the propagators for Tree-Level-Heavy (TLH)/Massless (TLM) 
modes split into a vacuum and thermal part as follows: 
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denotes the Bose-Einstein distribution function. 
TLM modes carry a color index 3 while TLH modes have a color index 1 
and 2. It should be noticed that the term cx u^u^, is due to the "propaga- 
tion" of the field, and = (1,0,0,0) represents the four-velocity of the 
heat bath [1]. 

The three- and four-vertices for the gauge vector bosons are represented 
in figure [2] and read: 
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Figure 2: Three- and four- vertices with Lorenz and color indices. 



The maximal off-shellness of momenta gauge modes due to the nontrivial 
ground state associated with the resolution |0| is constrained as: 



|p^-m2|<|0|2 (TLH modes) < |0|^ (TLM modes). (6) 



where for TLH modes the mass is given hy m — 2e|0| = 2ey ^ and for 
TLM wc have m = 0. 

The other kinematical constraint is on the center-of-mass energy flowing 
into a four- vertex that should not be greater than the compositeness scale |0| 
of the effective theory. For the momenta modes p and k entering the four- 
vertex one has 



(4) 



[4] abed 



-ie^{2'Kf5{p + q + s + r) 
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(5) 
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\{p+kr\<\<p\'. 



(7) 



This relation puts a strong restriction on loop integrations. Now con- 
sider (pi,p2) and {p3,P4,) as the pair of ingoing and outgoing momenta as 
represented in the diagram below. 




The Mandelstam variables are defined as: 

s = (p^+p^f = (p^+p^y^ t= {p^-pif = {pi-p2f, 

U = {pA-piY = {P3-P2f- (8) 

Accordingly, relation ([7]) reads then in terms of s-, t- and u-channels as 
the following 

|(Pi+P2)^| < |0P (s-channel), Kps - < |0p (t-channel), 
|(P3-P2)'| < (u-channel). (9) 

For the three-vertex conditions iQ are already contained in ([6]) by mo- 
mentum conservation in the vertex. 

From the above conditions one can immediately see the increase in the 
number of compositeness constraints by the s-, t- and u-channels. Conditions 
(EI) and ([9]) imply that the higher the loop order, the more suppressed their 
contribution to a thermodynamical quantity. General arguments suggest 
that, apart from diagrams associated with one-particle irreducible resumma- 
tions of propagators, only a finite number of diagrams contribute to the loop 
expansion The aim of the present article is to demonstrate the numerical 
evidence of this at three-loop level. 
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3 Constraints and Compactness 



The following concerns the most important poinlQ about irreducible three- 
loop diagrams. In [[Ij the only irreducible three-loop diagrams were consid- 
ered, 




a) b) c) 



Figure 3: Irreducible three-loop contributions to the pressure. Solid (dashed) 
lines are associated with the propagators of massive (massless) modes. 

and their moduli contribution to the thermodynamical pressure were 
computed. For all these irreducible three-loop diagrams we have k < k, 
where = 6 is the number of radial independent loop variables {p^, for 
i = l,2, 3, and k = 7 is the number of constraints on them, which is counted 
as follows. For diagram (a) and (b) there are 3 compositeness constraint^ 
emerging from the effective theory over s-, t- and u-channels at the four ver- 
tex \{pi+P2y\ < liPs—PiYl < |(P2— Ps)^! < |0p and 4 on-shellness 
relations = m^, P2 = m^, p^ = m? andp4 = {pi+P2~PzY — "^^i ^^e latter 
is due to momentum conservation at the four-vertex. Thus, for diagram (a) 
and (b) we have = 3 + 4 = 7. Diagram (c) is subject to case differentiation 
by momentum conservation at the four-vertexEI. For the case where both of 
the massless modes propagate off-shell there are 5 compositeness constraints 

\{Pl+P2f\ < \<P\\ \{P,-Pl?\ < \<P\\ \{P2-PZ?\ < \<P\\ \P\\ < \<P\\ \pI\ < \<P\\ 

and 2 on-shellness relations P3 = m^, p\ = (pi+P2~P3)^ = making k = 7. 
For the case where one of the massless modes propagate off-shell while the 
other one is on-shell there are 4 compositeness constraints +^2)^! < 

^this was first pointed out in [5] 

^these are used in the next section to verify tlic supports of radial loop integration - 
integrations over (r,9,ip). 

^hy using p| = {pi + P2 — Ps)^ , one can find dimensionless equations in terms of radial 
and angular variables which contradict p\ = p\ = ^ and agree with either = 0, p| 7^ 
or pI = Q,pI^Q 
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\{P3 -Pif\ < \{P2 -Psfl < \plor2\ < and 3 on-shellness re- 
lations pI = w? , pI = {pi + P2 — PsY = rn?, Pior2 — making again k = 7. 
Thus, for both cases in diagram (c) we have fc = 5 + 2 = 4 + 3 = 7, and 
therefore we conclude that for all the irreducible three-loop diagram k < k. 
The fact k < k shows that irreducible three-loop integrations are either com- 
pact or empty [T], and this is in sharp contrast with the two-loop case. For 
instance, consider the following nonvanishing two-loop diagrams. 




Figure 4: Non-local and local two-loop contributions to the pressure. 

For both of these two-loop diagrams one has k = 4 since in each case 
{p^, \p \)i for i = 1,2. For both of these diagrams there is one composite 
constraint +^2)^! ^ |0P and two on-shellness conditions pi = m? and 
p\ = making = 3. Obviously, this implies that two-loop integrations 
are not compact [2], since k > k. 

The fact that k < k for irreducible three-loop diagrams gives a very 
strong indication for the (rapid) convergence in the loop expansion. As the 
numerical results of the next section will show, irreducible three-loop integra- 
tions have much more suppressed contributions to the pressure than two-loop 
integrations for which k > k. 

4 Numerics and Results 

The following explains the statistical method (Monte-Carlo) in subsection 
14.11 which is used then with the constraints emerging from the effective the- 
ory (see equations (14), (17) in pj) to verify the supports for radial loop 
integrations in the following subsections 14.21 and 14. 3[ Subsection 14.41 tests the 
reliability of the statistical method and compares the shape of the region of 
radial loop integration for two-loop with the shape of the region of radial loop 
integration for irreducible three-loop and comes to an important conclusion 
precisely matching the general discussion in section [3l 
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4.1 Monte-Carlo-Integration 



The following section describes how the integral of a function / over a region 
G C M" can be calculated using statistical methods. The region G is deter- 
mined by a set of inequalities. Its characteristic function will be denoted by 
Xg- The integral which shall be calculated reads: 

/ f{x)dx^ / XG{x)f{x)dx (10) 

Jo JR^ 

If G is compact it can be included in a box B — [xi, Xi] x [2:2, X2] x ... x 
[xn,X„]. Therefore the integral can be written: 

/ XG{x)f{x)dx ^ / XG{x)f{x)dx (11) 

JR" J B 

If B has volume V, yXs can be considered as the probability density 
function of a random variable X, which is equally distributed on the box B. 
The integral is exactly the expected value E{VxG{X)f{X)): 



j^XG{x)f{x)dx = V j^XG{x)f{x)^XBdx = VE{xG{X)f{X)) (12) 

The Monte-Carlo-Method to calculate this integral consists of a statistical 
estimation of the expected value. It is known that the mean value X of a 
sample is an unbiased estimator for the expected value. If one draws a 
random sample (xi, X2..., Xn) of points from the box the estimation becomes: 

1 

f{x)dx^V-y^XG{xi)f{xi) (13) 

1=1 

The Monte-Carlo-Method is particularly useful to determine integrals 
over high-dimensional integration regions, where deterministic methods would 
be too time consuming. Unfortunately, the precision of this estimation in- 
creases only with the root of the sample size. This means that the Monte- 
Carlo-Method is a fast way to achieve a result with a relative precision of 
about 1 percent, but to achieve one more decimal place in the result the 
sample size must be increased by 100. 
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4.2 Computed ratio of |APa| and lAP^I to Pi^ioop 

The calculation of APa and AP^ in fj] includes an integration over a six- 
dimensional region G with variables of integratioE0 (xi, X2, X3, Z12, -Zia, -223), 



223, > 



dxi dx2 / (ixs / ^2:12 / dzis / ^2:23 x 

l,m,n=r J J J J Jz23,i 

f{xi, Zij, A, e) 

5(4e2 + e)5(x2, e) - (14) 

(-l)'+M(xi,e)C(a;3,e)- 

(-l)-+"i?(x2,e)C(x3,e) + ...), 



where A, B and C are square roots a/x? + 46^ and " ..." represents linear terms 
in Xj and Zij for z,j = 1,2,3 with A and e(A) as dimensionless temperature 
and effective gauge coupling, respectively. The whole integrand results from 
contracting the Feynman rules explained above. The Delta-function is left 
explicit in order to emphasise its non-trivial and precise integration later. 

As the Zij stand for cosine values, it is possible to restrict their range to 
the interval [—1,1]. The exact shape of G is determined by the (rescaled) 
compositeness constraints 

2:12 < -^{Ae'^ - A{xi,e)B{x2,e) + = g{xi,X2), 
X1X2 2 

^13 > ^(-4e2 + A(xi,e)C(x3,e)-^) =^(xi,X3), (15) 
X1X3 2 

^23 > -^{-Ae^ + B{x2,e)C{x2„e) -]-) = g{x2,x:i), 
X2X3 2 

and the additional angular condition on Z23,u{i) as defined in [1]. 

This angular condition allows a further reduction of the range of ^23 to 
the interval [z23,i, -223,u]- Obviously, this interval depends on the values of Z12 
and 213. 

The compositeness constraints and the additional angular condition allow 
a restriction of the Xj to the interval [0, 3] through considering the support 
for these radial independent loop integrations which does not contradict a 
nonempty support for the angular integrations [2] • Therefore a box B which 
contains G is identified. The compositeness constraints then define the char- 
acteristic function xg whereas the lower and upper limit for ^23 are used 
directly as integration limits. 



* Notice that xi is integrated analytically through the delta function in APa and APf,. 
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On the contrary to diagram a) and b) the region of radial loop integration 
for diagram c) cannot be determined in the same way. This has mainly to do 
with the complexity and increase in the number of compositeness constraints 
according to additional off-shell variables with arbitrary time components 
implicated by the four-vertex momentum conservation relation of diagram 
c) . Apparently, this makes it impossible to determine the absolute value sign 
of these compositeness constraints from which a small compact region for 
radial loop integrations could readily follow (respecting a definite nonempty 
support for angular integrations). The constraints on diagram c) are so re- 
strictive that its region of integration turns to be empty (as explained in 
the next section). According to the to the above analysis the ratio of the 
moduli APa and APf, to one-loop are depicted as a function of the dimen- 
sionless temperature Ac = 13.8 < A < 140 as the following. The one- loop 
pressure does not include the ground-state contribution. For the effective 
gauge coupling, the plateau value e = 8.89 is used for all temperatures A [3]. 
Throughout most of the deconfining phase, this is admissible, but the loga- 
rithmic pole of e at the critical temperature is ignored. Figure 0- El show that 
at the critical temperature Ac = 13.89 there is no loop contribution and AP^ 
and APb approach zero for large temperatures. This is due to the fact that 
with rising temperature the monopoles become massive and dilute, and the 
scattering processes are suppressed. At asymptotically high temperatures 
the contributions remain finite. The maximum of the ratio of the moduli to 
one-loop are peaked between A = 17.5 and A = 20 by 6 ■ 10"^^ and 2 ■ 10~^ for 
diagram a) and b), respectively. This shows the dominance of diagram b) in 
the irreducible three-loop expansion (the next section shows that diagram c) 
is vanishing). Comparing figure [6] with [7| and [HI shows the significant domi- 
nance of Coulomb fluctuationti (10^^) over quantum fluctuation (10 ^^) for 
the massless propagation in the contributions of diagram b). 



owe this particular analysis to discussions with Markus Schwarz, see [H [5] 
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Figure 5: Upper estimate for the modulus \ APa\/ Pi_iQQp as a function of A 
for diagram a). 
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Figure 6: Total upper estimate for the modulus | APbl/Pj^ j^^p as a function 
of A for diagram b). 
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Figure 7: Upper estimate for the modulus of the Coulomb part of 
loop ^ function of A for diagram b). 
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Figure 8: Upper estimate for the modulus of the quantum fluctuations of 
\APb\/P^_YQQ-p as a function of A for diagram b). 
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4.3 The (empty) region of integration for |APc| 

Consider a box with volume 10 x 10 x 10 x 20 x 2 x 2 x 2 for the loop and angular 
variables xi,X2,X3,yi, zu, zi^, Z23. Through condition lyf — xH < 1 in [1], and 
the definition of Z23,(u,i) a subset of this volume is determined, which roughly 
represents 2 % of the volume of the box. In this subset 150,000,000 points are 
chosen for Xi, X2, X3, yi, 212, -213, -^23 randomly and all four possible ± combi- 
nations for 1/2 and are taken into account. Then it is checked whether the 
conditions (1), (2), (3) and (5) are satisfied which leads to 600,000,000 tests. 
No point satisfies these conditions all together which estimates a fraction 
1/150,000,000 of the subset as the region of integration. This makes a fraction 
by 1/150,000,000 x 0.2 regarding the box. The box has a volume of 160,000 
which results a maximal volume of 1/150,000,000 x 0.2 x 160,000 ^2x 10"^ 
for the region of integration. The typical length distance is 0.2 which means 
v^2 X 10~^. It is then highly probable that the region of integration is 
empty. A similar analysis is applicable to the vacuum-thermal case show- 
ing that its region of integration is also empty. It is also possible to argue 
in terms of compositeness constraints that the region of integration for the 
vacuum-thermal case is automatically empty when the region of integration 
for vacuum- vacuum is empty [1]. 

4.4 Two-Loop Check 

The Monte-Carlo Method is now used to calculate the results for the local 
two-loop diagram b) in figurelH which has already been evaluated analytically 
[HI [7j . The results are compared in order to check the reliability of the Monte- 
Carlo Method. The calculation of APj, of figure 2 in [21 El El includes radial 
loop integration over a 3 dimensional region G with variables of integration 

(X, y, Zxy), 

j dx j dy J dzxyf{x,y,Zxy,e) x (16) 
nB(27rA-3/2^(x, e))nB{27(X-^/^B{y, e)) 

with 

z.y < — (4e2 - A{y, e)B{y, e) + ]-)= g{x, y), (17) 
xy z 

where the Bose factors are left explicit to emphasise their role in radial 
loop integrations as follows. 
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In order to carry out this integration with the Monte-Carlo method one 
needs to bound the region of radial loop integrations by restricting the Bose 
factors to the following interval : 

< 27rA-^/^A(x, e) < 10 (18) 

and 

<2nX-^/^B{y,e) <10. (19) 

These restrictions ensure that no Boltzmann tails associated with maxi- 
mal Boltzmann suppressions are included in the region of radial loop inte- 
gration for diagram b) in figure HI The above integrand is plotted in figure 
M below. 

200 




Figure 9: The integrand in ffTUj) is plotted as a function of x and y for fixed 
Zxy = COS Z{x,y). The horizontal plane represents the x-y-plane, where the 
integrand (mountain-shaped object) stands on the 'infinite' strip - the domain 
of definition for the contribution in figure fTOl . 
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In the x-y plane from x = y = 150 the maximal suppressions becomes 
evident which justifies the taken limits ffT^ and ffT^ . For instance, by taking 
x = y= 150, e = 8.89 and A = 30, the value of 27tX-^/'^A{x, e) is ~ 5.77. 

This shows that the restrictions in ( fT8l) and ( |T9l) are reliable estimations 
for bounding the region of radial loop integration. It should be noticed from 
figure - El and [TT] that the upper estimate for the contribution, which is of 
our main interest, takes its values for temperatures 20 < A < 30. 

There is an analogy between figure [9] and figure [10] in two dimensions. 

y 




1 2 3 4 5 6 j( 

Figure 10: Analytical determination of the integration limits for z^y = 
cosZ{x,y) depending on x and y. In the small triangular region near the 
origin, the upper limit of integration is 1; in the diagonal unbounded strip, 
the upper limit is g{x, y). From the two regions adjacent to the strip and the 
axes no contribution arises [7]. 

The interior region within the infinite strip of figure [H] is exactly the 
domain of definition of the integrand in the x-y plane of figure [TUl where 
from X = y = 150 the maximal suppressions become evident. 

The ratio of two-loop to one-loop computed by the Monte-Carlo method 
based on the above restrictions is represented in figure [TT] below. There are 
no significant statistical deviations about the analytical results computed 
in P, [7j. The upper estimate of the analytical and statistical results for 
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AP/Pi_ioop are both of the order 10^^ at temperature^ A = 20 ± 5. Ob- 
viously this striking compatibihty shows the rehabihty and correctness of 
the Monte-Carlo method for such radial loop integrations. It ensures that 
the irreducible three-loop results represented in this chapter are reliable and 
correct results. It should be noticed that these irreducible three-loop inte- 
grations are extremely sensitive and the slightest shift in the calculations 
can lead to significant changes in the numerical results. For instance, one 
should be aware of the simultaneous ± sign couplings of the algebraic prod- 
ucts within the arguments of the delta functions when summing over their 
zeros in integrations- not all free combinations contribute 



AP/Pi_ioop 



- 5 X 10"^ 

- 1 X 10"^ 
■ 1 . 5 X 10"*^ 

-2 X 10"^ 




20 40 60 80 100 120 140 



Figure 11: The computed ratio of APf, and Pi-ioop with the Monte-Carlo 
method plotted for 13.89 < A < 140. 



Let us now compare the integrand of (19) in [2] illustrated in figure M 
for the (local) two-loop diagram b) in figure H] with the integrand of (12) 
in [1] illustrated below in figure [I2] for the (irreducible) three-loop diagram 

^Notice that the plateau value of e in [H[7] is not 8.89 rather than 5.1. 

^For example, in products containing A, B and C terms one should not consider ex- 
pressions like d{ztAB ± AC ± BC) for which there are 2^^ = 8 free combinations. This 
is because any choice of sign for A, B and C in 'any one' of the products AB, AC and 
BC is simultaneously the same choice in the 'other ones' (simultaneous sign couplings). 
To be more specific, if the sign of A in AB is taken negative so it is taken in AC, there- 
fore we should consider expressions hke 5{{-l)''+''AB + (-1)°+'=^C + (-1)''+^BC) with 
a,b,c— 1,2 for A, B and C, respectively. The same sign coupling condition holds for 
such products within the unresolved compositeness constraint (17) which makes it so re- 
strictive. Taking all combinations free and ignoring the simultaneous sign couplings is 
including points which do not belong to the region of radial loop integration and by that 
the contributions can be increased significantly. Therefore, the precise Delta Function 
integration is extra explained in the appendix. 
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b) in figure El To get an idea of tlie precise sliape of the region of radial 
loop integration for the three- loop diagram b) in figure [3] illustrated in figure 
[T2] the following procedure to reduce the number of dimensions has been 
adopted: 

A grid of cosine values for 2:12, -Zis, and ziz with width 0.15 is used. This 
amounts to approximatively 2300 combinations of these values. For each of 
these combinations the integrand is calculated as a function of X2 and X3. 
Then all the diagrams of these functions are superposed to get the upper 
envelope of these diagrams. Figure [12] shows this superposition. Figure [12] 
can be interpreted as follows: The area where the superposition (maximum) 
vanishes, does not belong for any combination of the cosine values to the 
region of radial loop integration. So the Xi — 0:3 extension of this region is at 
most the area, where the plotted upper envelope does not vanish. It can be 
seen that figure [I2] supports the assertion in [2] that the region of radial loop 
integration is bounded and hence compact. The fact that the upper limits of 
integration for X2 and X3 have been set to 3 ensures under all circumstances 
that the region is fully covereclf]. 

Recall that Figure [H] shows the integrand for the two-loop diagram for a 
particular choice of the cosine value Zxy close to -1. The region of radial loop 
integration is obviously not compact, since the support of the integrand is 
an 'infinite' strip in the x-y plane of figure [10] Nonetheless the integrand is 
suppressed the further x and y are away from the origin. 

The remarkable contrast between the two-loop and the (irreducible) three- 
loop diagram is again, as it can be seen from figure [9] and [T2l that the region of 
radial loop integration in the latter case is compact. This was the analytical 
claim of 



^This is exactly the claim in [2] that the support in X2, 2^3 for the integral in x.yz 
is contained in the compact set {x2,a;3,< 3} while the support for the integration in 
2:12, zi3, Z23 naturally is contained in the set {-1 < zi2, Z13 < +1; ^23,; < -^23 < 223,«}- 
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Figure 12: The integrand in (fT^ is plotted as a function of X2 and 
X3 for a grid of cosine values Z12 = cos /(aTi, a^), 2:13 = cos Z(xi, 0^3) and 
2:23 = cosZ(ar2,x"3). 

5 Summary, Conclusions and Outlook 

In the present work it was shown that the ramifications due to the increase in 
the number of constraints on the loop momenta of the irreducible three-loop 
diagrams emerging from the effective theory imply either compact or empty 
supports for integrations and accordingly very suppressed contributions to 
the thermodynamical pressure. This was particularly due to the increase in 
the number of compositeness constraints by the s-, t- and u-channels. 

The series of extremely suppressed results 10"'', 10~^^ and confirm a 
rapid convergence in the loop expansion of SU(2) Yang-Mills thermodynam- 
ics. The fact that an irreducible three-loop diagram vanishes exactly, where 
the thermal contribution terminates completely, was according to the ex- 
tremely restrictive compositeness constraints that could not be resolved for 
radial loop integrations implying an empty support for integrations. 

Comparing the modulus of the dominant irreducible three-loop contri- 
bution with the smallest two-loop contribution reveals that they are nearly 
compatible, and comparing it with the dominant two-loop contribution shows 
a suppression by a factor 10"'*. Apparently, the dominant irreducible three- 
loop contribution was significantly dominated by Coulomb fluctuations over 
quantum fluctuations for massless propagations. The other nonvanishing ir- 
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reducible three-loop contribution is suppressed by a factor 10~^^ compared to 
the dominant two-loop contribution. The differences between two-loop and 
irreducible three-loop integrations were analysed according to the relation- 
ship between the number of independent radial-loop-momentum variables 
and the number of constraints on them. It was shown that the region of 
radial loop integration for irreducible three-loop integrations is, in contrast 
with the (noncompact) region of radial loop integration for two-loop, either 
compact or empty. This was also shown by the illustration of integrands 
according to the constraints on integrations. 

The small numerical results computed for the irreducible three-loop in- 
tegrations were in agreement with the general expectations regarding these 
diagrams. These numerical results were computed by the statistical Monte- 
Carlo method and this method was explained and also tested for two-loop 
integrations. The compatibility between the statistical and the former ana- 
lytical results for two-loop integrations ensured the reliability of the statistical 
method for the computed irreducible three-loop integrations. 

Finally, it could be very interesting to continue with this work by consid- 
ering higher loop diagrams, such as the ones represented below, which can 
be related to the irreducible three-loop diagrams. 

The challenge would be to see whether any of them survive the further 
increase in constraints, and hence whether any of them still have infinitesimal 
contributions at all. Also, it would important to find a rigorous mathemat- 
ical proof for the conjecture that the dominance of constraints over radial 
independent loop variables implies compact loop integrations. 
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A Delta Function Integration 



The integral over the products of delta functions in giving an estimation for 
the irreducible three-loop contribution of diagram a) and b) reads as follows: 



dxi J dx2 J dx3 J dyi J dy2 J dys J dzu J dzi3 J dz2z x 

5{yl - X? - Ae^)5{yl - x\- \^)b{yl - - 46^) x (20) 
5(4e^ + yvy^ - xxX2Zvi - {y\y?, - 3:1X32:13) - (^22/3 - x^X'^z^'i)) 

The integration will be done numerically after the Dirac delta-functions 
have been integrated away. The general formula for the integration of a 
delta-function is: 

/ dx 6if{x,p))g{x,p) = X; Jff^\''\, (21) 

n{p) denotes the number of zeros of the equation f{x,p) =0. This number 
depends in general on the parameter p. The different zeros are then denoted 
by Xi{p). To integrate the first delta-function the equation yf — xf — 4e^ = 
has to be solved for yi. 



y,,, = {-irJxl + Ae^ (22) 



201) then reads: 



dxi / dx2 J dxs J dyi J dy2 J dy^ / ^2:12 / dzis / ^2:23 
1 

X 



|2(-l)^v/^FT4i^| 
6{yl -xl- Ae^ )6{yl - x j - Ae') x (23) 

5(4e^ + {-lyJxj + Ae^y2 - X1X2212 - 



{{-ly^xl + 4e2?/3 - ^1X3^13) - (2/22/3 - a;2a;32;23) 
Proceeding in a similar way with the next two delta-function yields: 
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J dxi J dx2 J dx3 J dzi2 J dzis J dz23 

2 2 2 



yyy^ ' ^ x 

S{Ae^ + {-ly+^^xj + Ae^^xl + Ae^ - XiX2Zu - (24) 



{{-iy+\ xl + 4e\ xl + 4e2 - x^x^z.^) 



((-l)^+ya;i + 46^x1 + 4e2 - X2X3^23)) 

To remove the remaining delta-function the following equation has to be 
solved for xi. One should be aware that squaring an equation is a non- 
injective operation so that the algebraic solutions need not necessarily be 
solutions to the initial equation. 



4e' + + 4ey a;2 + 4e2 - 

2^1X2-^12 - ((-1)^+^x2+46^x1 + 462 - X^XsZ.s) - (25) 

((-l)^'+yxi + 4e2^xi + 4e2 - ^2X3^23) = 
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